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*eF ILE**1D**UVXPOWCJ 
UU UU VV VV XX XX PPPPPPPP 000000 
UU UU VV VV XX XX PPPPPPPP 000000 
UU UU VV VV XxX XX PP PP 00 00 
UU UU VV VV XX XX PP PP 00 00 
UU UU VV VV XX XX PP PP 00 00 
UU UU VV vv XX XX PP PP 00 00 
UU UU VV VV XX PPPPPPPP 00 00 
UU UU VV VV XX PPPPPPPP 00 00 
UU UU VV VV XX = XX PP 00 00 
UU Vv V XX XX PP 00 
UU UU VV VV Xx XX PP 00 
UU U vv vv XX XX PP 00 00 
UUUUUUUUUU Vv XX XX PP 000000 
UUUUUUUUUU vv KX XX PP 000000 
LL IIIIII SSSSSSSS 
LL III] SSSSSSSS 
LL I] SS 
LL I] SS 
LL I] SS 
LL I] SS 
LL I] SSSSSS 
LL I] SSSSSS 
LL I] S$ 
LL I] SS 
LL I] SS 
LL I] SS 
LLLLLLLLLL 11111 SSSSSSSS 
LLLLLLLLLL IIIT] SSSSSSSS 
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AAAOOOAOOAOO 
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eccccceoO COO 
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( pisreey ; Detailed Current Edit History 
(4) 183 BTSePOWCI - conrubaes ** INTEGER*4 


i 14 
INTEGER*4 power routine 


. ete te eRe Se Ge Ge Ge 
eeeeeeeeeeeeeeeeeneee 


"S-SEp=1986 17:09:07 EATHNTL. SRESUVKPOUCS MAR: 1 


“TALE UVXSPOWCJ - COnPLED #* INTEGER*4 power routine 
SIDENT /1-006/ File UVXPOWCJ.MAR Edit: SBL1006 


PITITIITITITITI LILI TTT LILL LLL LITLE Li Liiiiiiitiiiiiiiiiil 


COPYRIGHT (c) 1978, 1980, 1982, 1984 B 
DIGITAL EQUIPMENT CORPORATION, MAYNARD. MASSACHUSETTS. 
ALL RIGHTS RESERVED. 


THIS SOFTWARE is FURNISHED A aa A_ LICENSE AND MAY BE USED AND sor ine 
ONLY IN ACCORDANCE WITH TERMS OF SUCH Tt AND WITH THE 


OTHER PERSON. NO TITLE TO AND OWNERSHIP OF THE SOFTWARE IS HEREBY 
TRANSFERRED. 


THE INFORMATION IN THIS SOFTWARE IS SUBJECT TO CHANGE WITHOUT NOTICE 
eeePORATion NOT BE CONSTRUED AS A COMMITMENT BY DIGITAL EQUIPMENT 


DIGITAL ASSUMES NO RESPONSIBILITY FOR THE USE OR RELIABILITY OF ITS 
SOFTWARE ON EQUIPMENT WHICH IS NOT SUPPLIED BY DIGITAL. 


(RRR RBRARAARARALASAASASASARRSRASRASALELASLASARR SALAS ASARRRRARAARASR RSA AAAS SE SO 


FACILITY: Language support Library - user callable 
"ABSTRACT: 


COMPLEX base to INTEGER*4 power. 
Floating overflow can occur. 

Undefined Seppnens tes “he can occur if 
base = (0.,0.) and exp <=0 


VERSION: 0 


yk dh ta 
Jonathan M. Taylor, 29-jun-77: Version 0 


Edit history for version 0 


SOOOSCOCSCOSOOOOSOSOOOOSOSSOOSOSOSOOOSOOOSOOOOSOOOOOOOOOOOOoOO 
COOOOCQOoOCOoOoOooOooCoooooooooooooooooooooooooooooooo 
COOOCOCCOCCOOCOOOOOSOSOSOOOSOOOOSCOOSOOOOOOOOOOOOOSOOOoOOoOO 
fn De a ee ee) ee eee www wa lw lolol olololololololololeololeloleleololololololeolo ole) 
MEW 0 OONOAOU EWN $9 OD NAUES WN O OONOAUES WN OOONOUS Wh 


BPE FANNIN NINIAIDPONPONPOPIPINPONPI NN 2 SP OP Os Oo es a 


td Nata 


xe 
cal 
AO 


pLeR ** INT 
; Detai 


SOoOoQoQoooooooooooo 2- 
OOoOooooooooooooooo 
pole fale lel ole lal elo lololololololo) 
WR OOONOUSWR—OOOon 


PAEQAXA MMMM 


oooocoo 


ed Current Edt wth cory =$Ep=1986 99:28:07 EAtHNTe. sacsuvaPoucy.mar;1 °° (By 


~SBTTL HISTORY ; Detailed Current Edit History 


Ha ™ CJ 
- enenge H$SERROR to MTH HSS IGNAL JMT 5-OCT-77 
hange FORSFLAG_ JACKET to MTHSFLAG., JACKET. TNH 17-July-78 
- - Fix bug giving divide fault, cr wrong results for 
some negative pgyers. A4so clean up comments. 
SPR 20 -Oct-7 
1-001 - —— e version ou a to_1 ond 7TH. UNDE XP 
THSK_UNDEXP S 07=-DE 
1-008 ~ ada "to the PSECT sairastaves 98° 22-DEC-78 
- Declare externals. 
- Use general mode iqaressing. "SaL 30-Nov-1981 
- monover all es of floating instructions for Microvax. 
1-006 = Correct name to OTS$SPOWCJ. SBL 31-May-1983 


+. ptsters. fet Version 0 of 9 


K 
he INTEGER*4 power routine 16-SEP=- 


ee e@eee 


SoOoOOoCOCCoOoCOOCOoOoOOoOCOOoOoOCoOooOOCOoOOOooOoOoOoOSoSS }r 
pole lo lo lelolololeol ole ole ole oleolol elo olololol oe ol eolol ole — mm .s 


SOOOOOCOOCOCOOOOSOOOOOOOOOOoOSoOS 


OOOO 000000000009 00090909 0909 09 SII NIN SINS NINO OOOO 
NAUS WN (OOO NAME WN O OONAUE WN OOONOUM 


ooo 


“SEP- 
~SBTTL DECLARATIONS 


INCLUDE FILES: 
EXTERNAL SYMBOLS: 
-DSABL GBL 
eth MTHSS$SIGNAL 
-EXTRN OTSSDIVC 
eEXTRN MTHSK_UNDEXP 
MACROS: 
EQUATED SYMBOLS: 

base = 4 

exp = 12 


OWN STORAGE: 


PSECT DECLARATIONS: 


AX/VMS Macro V04-00 
MTHRTL.SRCJUVXPOWCJ MAR; 1 


; Math error routine 
3; COMPLEX division routine 


base input formal - by-value 
exponent input formal - by-value 


-PSECT _OTSSCODE PIC,SHR,LONG, EXE ,NOWRT 


; program section for OTS$ code 


| 


yur se out 
-006 


te 
a 


ek ek at ts = os = ts = = SS SS SS SS SS — SS 9 Sh — — Ss 2 2 9 ss SS Ss SS Ls ts a SS 2s SS 8s 2 SS 2s ss + os YS oOo 
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-SBTTL OTS$POWCJ = COMPLEX*8 ** INTEGER*4 


or 
= 
“nn 
SS 
ov 
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SOOSCCSOSOOCOOOOSOOOOOOSOSOOOSOOOSOOOOSOO OOOO O OOOO OOSOOOOOOOOOOOoOoOo om 
cm 


vu- 
ma 


p++ 
; FUNCTIONAL DESCRIPTION: 


COMPLEX result = COMPLEX base ** INTEGER*4 exponent 
The COMPLEX result is given by: 


base exponent result 

any >0 PRODUCT (base * 2**i) where 
i is each nonzero bit in 
exponent. 

(0., 0.) <=0 Undefined exponentiation. 

not (0., 0.) <0 PRODUCT (base * 2**i) where 
i is each non-zero bit in 
texponent:. 

not (0., 0.) =0 (1.0, 0.0) 


Floating overflow can occur. 
Undefined exponentiation occurs if base is 0 and 
exponent is 0 or negative. 
CALLING SEQUENCE: 
Power wfc.v = OTSSPOWCJ (base.rfc.v, exponent.rl.v) 
INPUT PARAMETERS: 
The base input parameter is standard FORTRAN COMPLEX. 
The exponent input parameter is a oleae lLongword integer. 
Both input parameters are CALL BY VALUE. 


IMPLICIT INPUTS: 
NONE 

OUTPUT PARAMETERS: 
NONE 

IMPLICIT OUTPUTS: 
NONE 


FUNCTION VALUE: 
COMPLEX base ** INTEGER*4 exponent 

SIDE EFFECTS: 
SIGNALS SS$_ARITH with floating overflow hardware code if 
floating overfl 


rflow. 
SIGNALS MTHS_UNDEXP (82 = ° UNDEFINED EXPONENTATION') if 
base is 0 and exponent is 0 or negative. 


NOME AR 9 OD NAME WIN | © OO NAME WIN SO OD NAME WI OOONAUE WN" OVOOnNOUSwnr- 32 


SoOCooooQoooooooQOoooooo oooooooooooooQooOoCooooQoQoOoOoOoOoOoOoOoOCoOoOoOooOoO «x 
PUPP PVPS B® Be EB BB EE ENN NAHI PO PONPINIPONPIPPOIYDDY 2 2 OO OOM MF OOOOOSOOCOCO O- 
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UVXSPOWCJ 16-SEP-1984 AX/VMS v04-00 
1-006 OTSSPOWCJ = COMPLEX*8 ** INTEGER*4 6-SEP=1984 2 94 93:97 :07 CMTHRTL. SREIU UVXPOWCJ MAR; 1 ° 3) 


| 
007C¢ 159 ENTRY OTSSPOWCJ, “M<R2,R3,R4, R5, R6> | 
160 ; “disable integer overflow 
2? 04 AC 7D 161 MOVQ base(AP), R2 ; R2/R3 gets COMPLEX base 
4 oO ag D 196 MOVL exp(AP), R4 : R4 = Longword exponent 
g 1 A 16 GEQ Ss 3 R4 = | exponent | | 
3% 4 cf C 164 MNEGL R4, RG | 
0c 54 98 § F 165 1$: BBCC #0, R4, EVEN 3; branc he if even and clear low bit 
50. 5 D 0013 196 MOVQ : RO/R1 = initial result 
54 54)COF FOF 6 0016 16 ROTL #-1, R4, RG ; R4 = unsigned. exponent / 2 
5B 1 01iB 168 BEQL DON : done if exponent was 1 
2A sii B20 198 BRB SQUAR1 3; else use rest of exponent 
RE 171 EVEN: 
QO1F fe 
OO1F 173; 
QOO1F 174 ; MOVD #1, RO ; RO/R1 = initial result 
Baie 13: 3; now becomes 
51 00 00 OO1F 177 MOVL #0, R1 3 Jnaginery part of initial result is 0 
50 08 50 b058 178 MOVF A 3; Real part of initial result is 1 
54 54 FF BF 9C 00 179 ROTL #-1, R4, RS ; R4 = unsigned_exponent / 2 | 
1D ig boSe 180 BNEQ SQUA : branch if exponent not 
§ 53 002C 181 TSTF R2 3 exponent was 0, text RP(base) 
4 i it} 86 BNEQ DONE ; done if non-0, answer is 1.0 
53 53 0030 18 TSTF R3 : IP(base) better not be zero 
44 12 ie 184 BNEQ DONE : it isn't return 1.0 
0034 185 
0034 186 UNDEFINED: | 
50 O01 OF 79 0034 13 ASHQ #15, #1, RO ; return RO = reserved operand 
7E 00° 8F 9A 0038 188 MOVZBL #MTHSK ONDEXP -(SP) ; FORTRAN error number 
00000000'GF 01 FB 44" 189 CALLS , G*ATHS$$SIGNAL 3 convert to 32-bit condition code 
004 190 ; and SIGNAL MTHS$_UNDEXP 
04 0043 191 RET 
0044 135 
0044 193 SQUAR: 
54 54 FF BF 78 +743 Ie ASHL #-1, R4, R4 ; R4 = treduced exponent: / 2 
0049 196 ; R2/R3 = square current base 
0049 197; 
0049 198 SQUARI: 
s$.6USSC2 B34? 199 MULF3 Re. R3, R5 ; RS = tmp = RP(base) * IP (base) 
26 2¢ 44 004D $30 MULF fk : R2 = RP(base)* 
5 5 44 O338 01 MULF ; R3 = iPibese}ee 
26 53 of B82 8 SUBF nS: ne : a = RP(base)**2 - IP(base)**2 
ss - 5 554 56 ADDF3 =R5, R5, R3 ; R35 = 2*RP(base)*IP(base) : 
E754 €9 ts Be BLBC R4; SQUAR ; branch if next exponent bit is 0 | 
33 88 ; RO/R1 = partial result * current power of base | 
55 53 50 45 0030 «308 MULF3 RO, R3, RS ; RS = tmp = RP(part) * IP(base) 
29 52.44 36) 09 MULF ‘ : RO = ng A RP (base) 
% 6 CS 51 45 64 10 MULF3 R11, R3, R6 3 R6 = = IP(part) * “TP (bas 
50 28 42 0068 11 UBF ; RO=RP (Bert) RP Cbase)-1P IP(part)*1P (base) 
51 § 44 Ree \¢ MULF R2, Ri ; R1 = IP(part) * RPC 
SS $0 of 1 ADDF R5, R1 ; RI=IP (parts sRP(base) sAP(part)#1P(base) 
54 54 C« RE CBE? 6 14 ASHL = #1, R4, RO : Ra = reduced exponent: / 2 | 
D1 = 12 76 15 BNEQ SQUAR1 ; loop if more exponent bits left 


N 14 
UVXSPOWCJ = COMPLEX ** INTEGER*®4 power routine = SEP-1984 mae Macro V04-00 Page 
1-006 OTSSPOuE = COMPLEX*8 ef INTEGER*4 -SEP-19 1382 oF: 96:33 YATURTE L.SRCJUVXPOWCJ MAR; 1 ° 
8 16 DONE: 
OC A Dd5 7 1 TSTL exp (AP) 3; test exponent sign 
1 78 18 BGEQ POWCJ ; done if positive 
5 5 7D 1 TSTF 3; test RP(result). 
2 1 7F 0 BNEQ RgCciP : if non-0, OK to take reciprocal 
a. 3 1 1 TSTF : RP(result) was 0, test [P(result) : 
AF 4 ; ; eacie BEQL UNDEF INED 3; undefined (0.0+0.01) ** <n 
m 3 FF 4 : MOVQ RO, -(SP) ; second arg pair is divisor 
0088 § : S*#1, (SP) ; push (1.,0.) on stack 
pees f ; becomes 
7E 00 00 0088 9 MOVL #0, -(SP) 3 push tnentoge part on stack 
7E 08 50 0088 330 MOVF  S*#1, =(SP) : Stack has (1.0 — 0.0) 
QOOO00000'GF 04 #«FB 43 H nents CALLS #4, G*OTSSDIVC ; RO/R1 = SRE 
04 0095 238 RET ; result in RO/R1 ; 
0096 34 
0096 235 . END 


a 
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Symbol table 


BAS = 4 
Bi ae 
EVEN iF R 1 
EXP = ¢ 
MTHS$S1GNAL tereerre x 
MTHSK_UNDEX eereetee 8X 
OTSSOTVC eteereeer § 6X 
TSSPOWC J RG 
CJ & 1 
RECIP 85 R 1 
saa ee 
UNDEF INED 800 32 8 1 
beer eon eon moran e = + 
H Psect synopsis ! 
PSECT name Allocation PSECT No. Attributes 
‘ ~*~ 00000000 < -} 00 ¢ 0.) NOPIC USR CON ABS LCL 
_OTSSCODE 00000096 ( 150.) O1¢ 1.) PIC USR CON REL LCL 
tow mn m nn wan wmeno nr neem ewan ane + 
: Performance indicators ! 
Phase Page faults CPU Time Elapsed Time 
Initialization 4 00:00:00.09 0:00:00.53 
zonnend processing 139 ah Sha 00:00:03.45 
ss : : ° : $ e 
Symbol table sort 0 0; 0:00:09 BR: 8889 «07 
Pass 2 56 Bat 0:00.46 ee? 64 
Symbol table output :00:00.0 00:00:00.05 
Psect synopsis output Ba! 338 °8 00:00:00.13 
Cross-reference output 8 O88 -0 BR: 88:00 08 
Assembler run totals 296 00:00:01.69 00:00:08.5 


pages) of virtual memory were used to buffer the intermediate code. 

pages of symbol table space allocated to hold 13 non-local and 1 local symbols. 
35 source Lines were read in Pass 1, producing 11 object records in Pass 2 

pages of virtual memory were used to define 0 macros. 


The work ing get Limit was 900 pages 
10 
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: Macro library statistics : 
Macro Library name Macros defined 
-$255$DUA28: CSYSLIBISTARLET.MLB;2 0 
0 GETS were required to define 0 macros. 


There were no errors, warnings or information messages. 


"SSEp=19be 14:29:07 


o V04-00 Page 


AX/VMS Macr 
SRCJUVXPOWCJ MAR; 1 


MTHRTL. 


NOSHR NOEXE NORD NOWRT NOVEC BYTE 
SHR EXE RD NOWRT NOVEC LONG 


— 


| 


WE Ecro nun statistics ~ (OULEK ** INTEGERSS pover rourhne "SSE TSRS 4:98:0% YATANTE. SRESGvAPOuey.mae:1 P29 8 


MACRO/ENABLE=SUPPRESSION/D1I SABLE=(GLOBAL , TRACEBACK) /LIS=LIS$:UVXPOWCJ/OBJ=O0BJ$:UVXPOWCJ MSRC$:UVXPOWCJ/UPDATE=(ENHS:UVXPOWCJ) 


| 
' 


| 
| 


| 0265 AH-BT13A-SE ) NT CORPORATION 
gave VAX/VMS V4.0 D PROPRIETARY 


